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This paper deals with the uniform exponential stabilities (UESs) of two hybrid control 
systems consisting of wave equation and a second-order ordinary differential equation. 
Linear feedback law and local viscosity, and nonlinear feedback law and interior anti- 
damping are considered, respectively. Firstly, the hybrid system is reduced to a first 
order port-Hamiltonian system with dynamical boundary conditions and the resulting 
systems are then discretized by average central-difference scheme. Secondly, the UES 
of the discrete system is obtained without prior knowledge on the exponential stability 
of continuous system. The frequency domain characterization of UES for a family of 
contractive semigroups and discrete multiplier method are utilized to verify main results, 
respectively. Finally, the convergence analysis of the numerical approximation scheme is 
performed by the Trotter-Kato Theorem. Most interestingly, the exponential stability 
of the continuous system is derived by the convergence of energy and UES and this 
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is a new idea to investigate the exponential stability of some complicate systems. The 
effectiveness of the numerical approximating scheme is verified by numerical simulation. 


Keywords: Hybrid system; nonlinear feedback; anti-damping; exponential stability; finite 
difference. 
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1. Introduction 


The wave equation on a segment with dynamical boundary control and local viscos- 
ity damping is a simple hybrid system, which is a mixed set of partial differential 
equation (PDE) and ordinary differential equation (ODE). More precisely, if let 
x € (0,1) be spacial variable, and t > 0 the time variable, the hybrid system is 
described by 


nuls, t) = Or. u(x, t) — a(x)d,u(a, t), (1.1) 
m0,,u(0, t) — 0,u(0,t) = F(t), (1.2) 


where, m > 0 is the tip mass, and for simplicity, the wave speed and the length 
of the segment are chosen to be unity. a(x) and F(t) are real scalar functions 
satisfying some properties which will be specified later and F(t) is the control at 
the boundary x = 0. It is well known that (1.1) and (1.2) can be used to describe 
both the vibration cable with a tip mass attached at the free end (e.g. see Ref. 21, 
Ref. 6, Ref. 20) and the mechanical oscillation inside the drill pipe (e.g. see Ref. 26, 
Ref. 25). 

In the case of a(x) = 0, designing suitable and effective feedback control laws 
to stabilize (1.1) and (1.2) were main contributions of the works mentioned above. 
For instance, Morgul considered the clamped boundary u(t, 1) = 0 at one extreme 
and linear feedback control law F(t) = —au,(0,t) + aust(0,t) with positive a and 
a at another extreme. The resulting closed-loop system is 


nulz, t) = Ona, t), 
u(t, 1) =0, 
mO;,u(0, t) — 0,u(0, t) = —auz(0,t) + auz(0, t). 


The exponential stability of (1.3)-(1.5) was lengthly discussed in Ref. 21. Moreover, 
Guo and Xu determined the optimal decay rate of the energy for a case left unsolved 
in Ref. 21, and the asymptotic expansion of the associated semigroup was also ob- 
tained ©. Recently, Vansprangh, Ferrante, and Prieur designed a class of nonlinear 
stabilizing feedback that only depends on the velocity at the controlled extreme. 
More precisely, the nonlinear feedback law F(t) = F(uz(0,t)) and nonlinear dis- 
sipative Neumann velocity feedback u,(1,t) = —g(uz(1,t)) were incorporated into 
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(1.1)-(1.2) and the corresponding closed-loop system was 


Onu(x,t) = srul, t), (1.6) 
mO,,u(0, t) — 0,u(0, t) = F(u:(0,t)), (1.7) 
Ux (1, t) = —g(ur(1,t)). (1.8) 


The exponential decay of the above system was investigated under different as- 
sumptions on the functions F and g in Ref. 25. The stabilization and regulation 
of a variant of hybrid system (1.1) and (1.2) was studied by using a proportional 
integral boundary controller in Ref. 26. However, the system of Ref. 26 is linear but 
the elasticity of the propagation medium is allowed to be non-homogenous. Further- 
more, if external disturbance was taken into consideration and entered the system 
through (1.2), Mei constructed an infinite-dimensional extended state observer to 
estimate the total disturbance and state simultaneously. An estimated state based 
controller was then designed to stabilize the hybrid system (1.1) and (1.2) 7°. 


Table 1. Known Results for Uniform Exponential Stability 


System Discrete scheme Verification method Y/N 
Wave Finite difference (FD) Spectral analysis®:24 N 
Wave Finite difference (FD) Numerical experiment? N 
Wave Finite element (FM) Spectral analysis? N 
Wave Finite element (FM) Numerical experiment? N 
Wave Mixed finite element Matrix inequality? Y 
Wave Mixed finite element Lyapunov function® Y 
Wave FD with viscosity discrete multiplier?’?4 Y 
wave FD with order reduction Lyapunov function!®30 Y 
Thermoelasticity FM Frequency domain!® Y 
Heat-wave coupling FD with order reduction Lyapunov function?’ Y 
Schrodinger on H'(0,1) FD with order reduction Lyapunov function! Y 
Schrodinger on L?(0,1) FD with order reduction Frequency domain” Y 
Eulur beam Finite volume Lyapunov function Y 
Timoshenko beam FD with order reduction Frequency domain?” Y 


In this paper, we continue to investigate the hybrid system (1.1) and (1.2) from 
the viewpoint of numerical approximation. That is to say we first need to perform 
numerical discretization for the component of PDE and transform it into a family 
of ODEs. A family of semi-discretized systems, which consists of pure ODEs, is 
derived from hybrid system mentioned above. We want to study the UES of the 
discrete systems, i.e., whether or not the discrete energies uniformly exponentially 
decay with respect to the mesh size. To draw a conclusion for this problem, one 
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should face three challenges. The first one is which approximating scheme is chosen 
among many approximating methods to discretize the spatial variable. The second 
one is which method is used to verify the UES. The last one is how to verify the 
solution of semi-discretized system converges to the solution of the original contin- 
uous system. Though one may can’t obtain uniform decay by using the classical 
finite difference and finite element schemes since these semi-discretizations produce 
spurious oscillations of the high frequencies 743!. The first question and the third 
question are easy to cope with because there are rich results on numerical solution 
to PDE !”. The second challenge is the most difficult one. However, some successful 
pioneering works have been given in last two decades. We refer to these references 
= and collect some related results in the table 1. 

The natural idea of verifying the UES is to follow each step of continuous case 
(e.g. see Ref. 5, Ref. 13-Ref. 15 and Ref. 24-Ref. 25). Whereas the method for the ex- 
ponential stability of the continuous may not suit for the discrete case or vice versa. 
Therefore, it is necessary to obtain UES without prior knowledge on how to verify 
the exponential stability of continuous system. On the contrary, the exponential 
stability of the original continuous system is obtained by UES of discrete systems 
and convergence analysis (see Theorem 2.4 and Theorem 2.5). This not only can 
greatly enlarge the study scope of UES but also conversely research the exponen- 
tial stability of the continuous system. For instance, to the best of our knowledge, 
the famous Huang-Priiss’s characteristic condition of exponential stability (e.g. see 
Ref. 16, Ref. 1, and Section 8.1 of Ref. 11) is difficult to be appled for the hybrid 
system with linear feedback law discussed in this note. But it’s discrete version is 
effective in the UES of semi-discretization systems (see the proof of Theorem 3.1 
below). In that case, we investigate the hybrid system with nonlinear feedback law 


and anti-damping by similar idea. Because the absence of characteristic condition 
a mman of exponential stability for nonlinear semigroup, we discuss UES by using the dis- 
£ = crete multiplier method. Both boundary and in-domain anti-damping, which are 
ad more interesting in terms of applications and were widely discussed,3:22-4:8-19,28,26 
are considered in our model. 
Thus, the contribution of this work is three-fold: 


e Extend the result of uniform exponential stability from PDE with linear 
feedback and boundary damping to the PDE with nonlinear feedback and 
anti-damping. 

e The characteristic condition of uniform exponential stability for a family 
of semigroups of contraction is successfully applied to the discrete system 
of boundary control hybrid system. This method has potential applications 
for the system in which the Lyapunov function is absent. 

e Present a new idea to discuss the exponential stability of distributed pa- 
rameter system, i.e., the uniform exponential stability and convergence of 
the energy implies that the continuous system is exponential stable (see 
Remark 2.1). 
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The structure of this paper is organized as follows. In section 2, we study the 
hybrid system consisting of (1.1) and (1.4) with linear feedback (1.5). A suitable 
semi-discretization system, which is uniformly exponentially stable with respect to 
the discretized parameter and convergent to the continuous system as mesh size 
tends to zero, is built. The exponential stability of the continuous hybrid system 
are also discussed in this section. In section 3, similar results of the second hybrid 
system consisting of (1.1) (1.7), and (1.8) are obtained. The effectiveness of the 
numerical approximating algorithms are also verified by numerical simulations in 
this section. In section 4, concluding remarks and future directions are shown. 


hag 2. The hybrid system with linear feedback 


In this section, we present main results on the hybrid system (1.1), (1.4), and (1.5). 
Firstly, we transform them into a first order port-Hamiltonian system with dynami- 
cal boundary condition, which is a natural treatment for the wave equation (e.g. see 
Ref. 19 and Section 7.1 of Ref. 11), and discretize the resulting system by average 
central-difference method. Secondly, the UES of the discrete systems is proved by 
the frequency domain method. Thirdly, the stability and the consistence of numer- 
ical approximating scheme are presented by using the Trotter-Kato Theorem. In 
the sequel, we denote A* in stead of A! the transpose of a matrix A. The prime ’ 
denotes a differential with respect to the temporal variable. Moreover, let N € N 
be a positive integer and h = 1/(N + 1) step size, the index k and j are fixed to 
denote k = 0,1,---,N+1 and j =0,1,--- , N throughout the whole paper. In this 
section, a(x) is a bounded nonnegative continuous function on (0, 1), i.e., a(x) > 0 
for x € w and a(x) = 0 for x € (0,1) \w with w C (0,1) being a nonempty open set. 


— 2.1. Order reduction and discretization 


Firstly, we define the following intermediate variables u(x,t) = Oju(z,t), w(a,t) = 
O,u(z,t). Then, the hybrid system consisting of (1.1), (1.4), and (1.5) can be rewrit- 


ten as 
v(x, t) = 0,w(a,t) — a(x)v(t, x), (2.1) 
Opw(a,t) = O,v(2,t), (2.2) 
n' (t) = w(0, t) — av(0,t), (2.3) 
v(1,t) = 0, (2.4) 
v(x,0) = v(x), w(x,0) = w(x), (2.5) 


in which 7(t) = mv(0,t) — aw(0,t), v°(x) and w°(x) are the initial velocity and 
the initial force of the elastic cable. Define the state space of (2.1)-(2.5) to be 
X = L?((0,1);C?) x C with the inner product 


~ pas 1 = 
(Y, Y)x = (v, ¥) r2 + (w, W) r2 + ———nh 
m+ aa 


January 16, 2023 8:23 WSPC/INSTRUCTION FILE ws-m3as 


6 F. Zheng, L. Zhang, Y. He 


for Y = (v,w,n), Y = (č, 0,Ñ) € X and (.,-);2 is canonical inner product on 
L? (0,1). Now we can formulate (2.1)-(2.5) into the abstract Cauchy problem 


Y'(t) = AY (t), Y(0) = (v°,w°,m)* € X, (2.6) 
here, Y(t) = (v(x, t), w(x, t), n(£))*, no = (0), and the operator A is defined by 


W(t, x) — a(x)u(t, x) 
AY(t) = Ve (t, £) (2.7) 
w(0, t) — av(0, t) 


and 
= D(A) = {(v,w,7)* € H*((0,1);C?) x C : v(1) = 0, n = mv(0) — aw(0)}. 
s (2.8) 
Using Lumer-Phillips Theorem, we can obtain the well-posedness result of (2.6). 
The proof is omitted since it is similar to that of Theorem 1 in Ref. 21. 
Theorem 2.1. The operator A: D(A) — X defined by (2.7) and (2.8) generates 
a Co semigroup of contraction T(t) on the state space. 
In fact, it is easy to see that for any Y € D(A), we have 
2Re(AY, Y)x = (AY, Y)x + (Y, AY)x 


0 w 


m + aa 


which implies that 


marroa- OU AE 


— m + aa 


/ a(ax)|v(a)|?da. (2.9) 


Fr However, it is well known that (2.9) has tight relation with the energy functional 
Y E(t) of (2.1)-(2.5). In other words, we define E(t) of (2.1)-(2.5) by 


E(t) = F) |v(t, x)|? + |w(t, x)| dz + ne (2.10) 


and differentiate it along the solutions of (2.1)-(2.5) to obtain 
E'(t) = Re(AY (t), Y(t)) x 


= - | a(z))v(0,)Pae - 


2 2 
m q aa elel, £)| + a|w(0, t)|*]. (2.11) 
Remark 2.1. Up to now, we don’t know whether or not the system (2.1)-(2.5) is 
exponentially stable. It is not easy to utilize the Lyapunov function methods of the 
Ref. 5, Ref. 30, and Ref. 21 and frequency domain method of the Ref. 1, Ref. 16, and 
Ref. 17 to verify this property. We could discuss the UES of (2.1)-(2.5) without prior 
knowledge of the exponential stability of continuous system. However, the UES and 
convergence of the energies between the discrete system and continuous implies that 
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the continuous system is exponential stable. This is one of main contributions this 
paper. 
Secondly, we introduce the so called average central-difference semi discretization 


scheme for (2.1)-(2.5). To this end, insert N + 2 points and N + 1 points, denoted 
by £k = kh and yj = (j + 1/2)h respectively, in the domain [0, 1]. The notations 


fii- fj fiait fi 
Oxfj4i = —— ae i= = 


denote the central-difference operator of f(x) and the average operator of f(x) at 
the node yj, respectively. If x in (2.1) and (2.2) is replaced by y;, then we have 
Ovv(y;, t) = Onw(y;; t) a aly; )u(y;, t), 
brw(yj, t) = OxU(Y;; t). 
In the above identities, approximating the temporal derivatives of the left-hand 


side by the average operators and the spacial derivatives of right-hand side the 
central-difference operators respectively, one derives 


Viya t) = dw; 41(t) — ajv,;41 (4), 2.12) 
Wit t) = 620; 41(¢), 2.13) 
nf (t) = wo (t) — avo(t), 2.14) 
n(t) = muo(t) — awo(t), 2.15) 
un41(t) = 0, 2.16) 
vk(0) = vg, we (0) = wk, 2.17) 


where v;,(t) and w;(t) are grid functions at grids z, and v? and wÌ the approxima- 
tions of the initial values v? (£) and w? (£p), respectively. Notice that the viscosity 
terms a(y;)v(y;,t) are in fact approximated by a(x;)[v(£j+1, t) + v(a;,t)]/2. Please 
reader to Ref. 5, Ref. 13, and Ref. 30 for more information. 

Take the state spaces of (2.12)-(2.17) to be 


X, = CNT x CNH x C. 


Assume that Y, = (v,w,n) E€ Xp with v = (vo,v1, vy) and w = 
(wi, W2,--+ ,wn41). If setting wo = a (mvo — 7) in light of (2.15), then we can 


define the inner product of X, by 


N = 
(Yn Fad = DD Pao aeae a e ae ay 


for Yp = (v, w,n), Y, = (0, 0,7) © Xa. 
Let Yn (t) = (up (t), wa(t), n(t))* with 
vn (t) = (vo(t), vi (t), , vn (2), (2.18) 
wht) = (wi(t), wa(t),--- ,wn4i(t)), (2.19) 
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being unknown variables of (2.12)-(2.17). The energy functional, which is denoted 
by E(t), of (2.12)-(2.17) is defined by 


v Ol + 3] + i (2.20) 


N 
En(t) = ; (Yn (t), Yat), = >> | 
j=0 


Finally, Æp (t) is non-increasing and this is parallel with (2.11). More precisely, we 
have 
Proposition 2.1. The discrete energy functional Ep (t) satisfies 


2 
E,(t)=—h a | 


jhew 


[ma|vo(t)|? + a|wo(t)|]. (2.21) 


tea om + aa 
T= Proof. Multiplying (2.12), (2.13), and (2.14) by hvj4ı (t), hwj4ı (t), and n(t)/ (m+ 
Q) aq) respectively and summing for j from 0 to N, we derive 
2 
= O=- D> a; lOl 
= jhEw 
n(t)n'(t) 


N 
Por p3 [rie Oswa (t) + wya (Hrv E). (2-22) 
= 


However, by some simple calculations, we have 


N 
>> [e+ rwg) 4 wiy (5243 (0) 
j=0 


a Uj+i(t) + vilt) (wj+1lt) — w(t y wy4i(t) + wi (t))(vj4i(t) — vilt 
D Orn oy +1(t) ONC (t)) 
=un+1(t)wn+1(t) — vo(t)wo(t). (2.23) 


Thus, (2.14)-(2.16) and (2.22)-(2.23) implies that (2.21) holds. 


2.2. Uniform exponential stability of (2.12)-(2.17) 


To present main result of this subsection, we first rewrite (2.12)-(2.17) as state space 
formulation. For this purpose, let e; = (1,0,--- ,0) € Ct, and 
11 1 
TN -1 1 
BNn+ı = , , Mn+ = S , 
—11 


Anyi = diag{1,0,--- ,0}, Ag = diag{ag,a1,--+ ,ay} belong to CN+YxN+), Set 


Bird 0 0 
Bh = 2 2ma~tAn41 By 41 —2a7 leš 
0 0 2 
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and 


—ma-'Ana1 = 271hAaBN41 MyN4+1 aleš 
Mi => Ves 0 0 
h(a~'m — aje 0 —ha™! 


Because B, is evidently invertible, (2.12)-(2.17) are then equivalent to the differen- 
tial equation 


Bry; (t) = MnYn (t) S Y; (t) = AnYn (t), Yn (0) € Xa, with An = B; Mh. (2.24) 


It is easy to see that B, is corresponding to the average operator of time deriva- 
tive of (2.12)-(2.13). If one replace Bp by identity operator, then the classical finite 
difference scheme of (2.1)-(2.5) is easily restored from (2.24), i.e., 


Yr (t) = MhYn (t), Yn (0) E X}. 


Here we explain the significance of the discrete scheme (2.24). We plot two figures 
in Figure 1- Figure 2, respectively. Figure 1 depicts the maximal real parts of Ap 
and Mp for N = 40 : 5 : 400. Figure 2 depicts the distributions of the eigenvalues 
of A, and Mp in which N = 500. 

We see that the real parts of the eigenvalues of Mp approach to zero and those 
of An approach to a negative number from these figures. In these figures, we take 
m=a=1, a= 2 and a(x) = cos((nx)/2). Numerical simulation results show that 
the classical finite difference scheme of (2.1)-(2.5) is not uniformly exponentially 
stable. This conclusion is consistent with earlier research results of Ref. 24. However, 
Figure 1 manifests (2.24) is perhaps uniformly exponentially stable. In the remaining 
part of this subsection, we give strict proof about uniform exponential stability of 
(2.24). 


x with the average operator 
0.2 F x without the average operator | | 
n 
o 
3 x x * Xx & X * XK XX KX * RK RN KR 
@ OF 2K OOK OOK OOK OK OOK OOK OOK OOK IIK OOK IDK IOK OOK OOK OOK OOK HOC DK 
S 
2 
D 
© -0.2 
2 
= 
g 
a 
T 
204+ 
T 
£ 
x 
© 
= 06 
-0.8 : ; 
0 50 100 150 200 250 300 350 400 


Fig. 1. Maximal real parts of eigenvalues 
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On the one hand, the inner product on the state space X, can be rewritten as 
(Yu, Yi) =h (BrYn. Bnn) 
Xa 


for Yp = (v,w,n), Yn = (©, 0,7) € Xp. Here (-,-) are canonical inner product of 
C2N+3 and wo := at (mvo — n) is applied in the last identity. On the other hand, 
Yh € (0,1), the operator An, which is discrete counterpart of A, are dissipative on 
the space X}. Because the dissipativity of Ap is key in the next part, we put it in 
the following property. 


Proposition 2.2. For all0 < h < 1, the operator Ap is dissipative and generates 
= Co-semigroup of contractions Th on the space Xp. 


Proof. For any Yp = (v,w,7) € Xp, we consider 


(AnYn, Yn)x, = (MnYn, BrYn) 


2 N 
=—h J aj vy; +h% [Paphos 43 + Drpönjt + 
jhew j=0 


and 


(Yn, AnYn)x, = (BaYn, MnYn) 


=-h y aj vj; 


jhEw 


2 Arp tm ay oT 
+A [vjr +0544 5e% 43 | + . 
j=0 


m + aa 


By using 7 = mvo — awo and some simple operations, we obtain 


2 ma|vo|? + alwol? 


l Re(AnYn; Yn)x,, =—h 5 Qj loja 
— jhew 


2.25 
m+aa i ( ) 


3000 


x without the average operator 
x with the average operator 


2000 f 


=} 
=} 
3 


=] 
S 
Ss 


imag parts of eigenvalues 
o 


-2000 F 


-3000 i 1 1 1 1 
3 -2.5 2 -1.5 ai -0.5 0 
real parts of eigenvalues 


Fig. 2. Eigenvalue distributions with N = 500 
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Moreover, for any 0 < h < 1, the spectral set o( Ap) of Ap is contained in the 
open left half-plane of C. This is main content of the following property. 


Proposition 2.3. For all 0 < h < 1, iR is contained in the resolvent set p(An) of 
An- 


Proof. If there exist G € R and nonzero Yp € Xa such that i8Y;, = AnYn, then it 
follows from (2.25) that 


i j 2 mal|vo|? + alwol? 
S 0 = Re (iBYn, Yn)x, =Re(AnYn, Ya)x, = —h 5y aj lezal [vo] [wo] l 


f m + aa 
jhew 


This means that h >? j,cu A; |Vj42 |? = 0 and vo = wo = 7 = O since 7 = mvo — awo. 
Noting that iGY = Ap,Y is equivalent to 

iBvj44 = OnWj42, (2.26) 
iBwy42 = OxVj41- (2.27) 


Setting 7 = 0 in the above identities, we have 
ibhv, — 2w, = 0 
—2v; +ibhw, = 0, 
It follows from Cramer’s rule that vı = wı = 0 since the determinant 


iBh —2 
—2 ibh 


|=- -4 


is nonzero. Similarly, we can repeat this process to prove v;41 = Wj+1 = 0 by solving 
(2.26)-(2.27) and using v; = w; = 0. This yields Y, = 0, which is a contradiction 
and we complete the proof. 


We can now state the following uniform stability standard, which is given in 
Ref. 1 or Ref. 17 and will be applied in the proof of Theorem 2.3. 


Theorem 2.2. Let h* > 0 and (Sn(t)) be a family of semigroups of contraction on 
the Hilbert space (Xn), and let (Ap) be the corresponding infinitesimal generators. 
The family (S),(t)) is uniformly exponentially stable if and only if the two following 
conditions are satisfied: 

(i) For all h € (0,h*), iR is contained in the resolvent set p(An) of (An). 

(ii) suPre(o,n*) ger IBI — An)" lee, <% 


We are in a position to give key result of this section. 


Theorem 2.3. The family (Tn(t)) is uniformly exponentially stable, that is to say 
there exist two constants M > 0 and w > 0 (both independent of h € (0,1)) such 
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that 

|T.(t)||n¢x,) < Me™, vt > 0. (2.28) 
Proof. The proof is based on Theorem 2.2. Notice first that, for all h € (0,1), 
the family (Th (t)) form a family of contraction semigroups (see Proposition 2.2). 
The fact that the family (Ap) satisfies condition (i) follows from Proposition 2.3. 
In order to show that the family (Ap) satisfies condition (ii) we use a contradiction 


argument. If the condition (ii) is false, then there exist 8, € R, hn € (0,1), and 
Zh, E Xn, such that 


l@BnInn — Ann) ZF IIx, 2 MZ, IIx, 


in which J), is the identity operator on Xp, n = 1,2,--+-, and the norm ||- ||x, is 
induced by the inner product of Xp. Let (i8nJn, — Ahn) Zk, = FẸ, and we have 


Fre ln, Z Aban = Arn) Fh Xan 
Setting YP = FY /|Fr lx; we obtain |Y} |Ix, = 1 and 


|U" |x, < nt, with U” = (ibnIn, — Ann Yp. (2.29) 
More precisely, for U? = ((",6",0") with C” = (GR, GP, CR) € CMF, €” = 
(ER, ER, ER 41) € C%*1, and 6” € C, we have 
N, 
h R 2 \a” |? 7 
n||2 — am n n 2 
IU” Ilx,,, = 5) 3 en + Sji | + NEET Int. (2.30) 


with € = a~! (mC? — 0”). It follows from Cauchy-Schwartz inequality, (2.29), and 
(2.25) that 


Re(U", Yn )x,, = Re (Ann Yhn Yn), 


= hin 5 Qj 


jhEw 


2 n|2 n|2 
ma|ve | T aļwg| < nt. (2.31) 


yn 1 
I+ m+aa 


Therefore, by using inequalities (2.30) and (2.31), we have 


feats 2 = O(n’), j=0,1, Np (2.32) 

Eal? = O(n), j= 0,1, Nns (2.33) 

5 a; lors = 0O(n™t), (2.34) 
jhew 

err =Oln-"), (2.35) 

|we |? = O(n"). (2.36) 


The notation |w?|? = O(n—') in (2.36) means that there is a positive constant C 
(independent of n) such that |w8|? < Cn~!. The other notation O(-) in (2.32)-(2.35) 
has the same meaning. 
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Furthermore, combing Bh„U” with U” = (ibnIn, — Ahn )Yp , we obtain 


iBntzya — Jn W541 + av; 1 = E (2.37) 
By rearranging (2.37) and (2.38), we derive 
iPr hrf ya — wy ya = —2w} +p}, (2.39) 
in which 
pt = hn [ajo Ha- (2.41) 


Set 7 = 0 in (2.39) and (2.40) and solve the resulting linear equation to arrive at 


vots) =! (ibam 2 — 20} + p 
wori) 4+ BRR 2 ibnhn) \— 205 + Anori)’ 


and which implies 
opal =OR™), [w541 =O(n™), (2.42) 


by (2.32)-(2.36), (2.41) and some basic calculations. Now (2.35)-(2.36) and (2.42) 
together yields 


2 

jor? = [2v8 ; -v| = O(n), 
2 

jwp? = 2w; — we] = O(n"), 


By letting j = 1 in (2.39) and (2.40), we can obtain 
W, P =O), wht, = O(a), 


juz]? = O(n), fw? = O(n"), 


by using the same method as above. Graduall 


y, we can derive that 


Weapon wh, =O), 


hold for all j = 0,1,--- , Nn and even more 


i Í = O(n“) (2.43) 


since h,(N, +1) = 1. It is easy to see from (2.35)-(2.36) that 
In” |? = mv? — awe |? = O(n"). (2.44) 


Finally, (2.43) and (2.44) imply that |Y? ||x,,, = O(n7') for all positive integers 
n, which is a contradiction to ||Y;? |x, = 1. The proof of the Theorem 2.3 is 


complete. 
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2.3. Convergence analysis 


The goal of this subsection is to show the approximating generators Ap, on the 
spaces X}, such that the Co-semigroups T},(-) generated by An, approximate T(-) 
in some sense. We will need the following assumptions: for every h € (0,1) there 
exist bounded linear operators P, : X > X, and Enr : X, > X satisfying 

(A1) There exist two positive constants Mı and Mo such that || En ||x < Mı and 
lPrllx < Mə, 

(Az) |En PhY — Y ||x > 0 as h > 0 for all Y € X, 

(43) PrEn = In. 


™ Theorem 2.4. /Trotter-Kato'®] Assume that (A1) and (A3) are satisfied. Then the 
= following statements are equivalent: 
(a) There exists a Ay E€ p(A) O NFL p(An) such that, for all Y € X, 


|En(Aoln — An) t PRY — (Aol — A)7'Y||x > 0, as h > 0. 
(b) For every Y € X and t > 0, 
|EnrTh(t)PhY — T(t)Y |x > 0, as h —> 0 
uniformly on bounded t-intervals. 


However, in order to prove Theorem 2.4 one faces several major difficulties. 
The most difficult one is the direct verification of the consistency property (a). It 
involves computation of the resolvents (AgZ;, — An)~', which in general is almost 
impossible. Therefore, the following property is useful. With the aide of it one can 
replace (a) by a condition involving convergence of the operators Ap to A in some 
sense of Ref. 10. 


ry Proposition 2.4. Let the assumptions of Theorem 2.4 be satisfied. Then statement 
= (a) of Theorem 2.4 is equivalent to (A2) and the following two statements: 

(C1) There exists a subset D C D(A) such that D =X and (AgI — A)“ D = X 
for aro >w. 

(C2) For all Y € D there exists a sequence (Y) with Yp € D(Ap) such that 


lim EnyY n = Y, lim Enr AhY n = AY. (2.45) 
h-0 h->0 
Proof of Theorem 2.4 Let xs be the characteristic function of a set S' and define 
the matrix F},(x) consisting of characteristic functions y(2,,¢;,,) (J = 9,1,--- , N) 
by 
X(ao,#1]> X(z1,£2]> PUR X(an an 41] 
F(z) = X(zo,x1]> X(z1,£2]> ` p hi X(zN,£N+1] 


X(ao,¢1)> X(z1,z2]3 (°° > X(T N,£ N41] 
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For any Yp = (v, w, n) € Xp with v = (vo, v1,- +- , Un) and w = (w1, w2, , WN+1), 
set wo = @7! (mvo — n) and define the extension operators Ep, : Xp > X: 
N 
dy o; 42) X(a5, xj+1] 
EnYn = diag(Fh (x), Fa(2), 1)BhYn = | Dj=0(Wj+3)X æ] | - (2.46) 
n 


Choose the dense subset D ê D(A)N(C?[0,1])? x C of X. For any Y = 
(v(-), w(), n) € D, set U = (v(zo), v(z1), +++, v(en))*, W= (w(21), w(zo), =, 
w(rn41))*, and Yn = (0,W,n)*. By using v(1) = v(xn41) = 0 and w(x) = 
a~'[mv(x9) — n] (it comes from the definition of the operator A), it is easy to see 


that 
Nf v(ej4i)+v(25) 
2-0 (—— | X[z;;£;+1] 
Y, = N (wlæj+i)+w(z;) 
EnYn par, (seer) Xizi zj] |? 
n 
and 


N w(x —w(x 
7 A ( eta): ( 2) Wee 
En AnY n = diag(Fh (x), Fh (2), 1)Mh Yn = | 5Y (Hee) Xepe 


{= 


mv(zo) — aw(zo) 


By using the same method of Section 4 in Ref. 10, we know that the conditions 
(C1) and (C2) of Proposition 2.4 hold. 

To prove (b) of Theorem 2.4 and give the main result of this subsection, we need 
further to verify (Aı)-(A3). In view of (A3) and (2.46), construct the projecting 
operators Ph : X > Xh 


©) Il fi(a) 


T°vu(a) T fil 
(x) 
PAY = | wla) |, PA) =2|) , 
n 
Infi) 
for Y € Xn, l = 0,1, fo(x) = v(a), fı(x) = w(x). Iu(x) and Ij ,w(x) are defined 
by 
Bo eu (x)dz, 
Io( =l Pal a)dx — I}, v(x) 
ae =f SE wl w(a)dz — ag *[I$u(ax) — n], 
Ty, w(x) nef a a)dx — I;w(x). 


In light of the definitions of E» and Ph, it is easy to see that they satisfy (A1) and 
(A3). To prove that Eyn and Py satisfy (A2), we firstly assume Y = (v(x), w(x), n) € 
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D(A). With this assumption we have 


N Tj 
2-0 Tj ! v(x)d2X]e; wjp] 
Ep P Y = h`! y j ae w(x)dzrX{e,;, 
hn 


x41] 


and 


FD) E v(T)] Xle; e41] 
E,P,Y-Y=| 5% [v(9%) — w]X{e;,0;41) | > 0 as h> o, 
0 


here 07 and 0% are chosen such that i v(a)da = v(0}7)h and | w(x)de = 
w(0})h when the mean value theorem is applied, and the continuity of v(x) and 
w(x) on the internal [0, 1] is applied in the last step. Thus combing this result and 
the density of D(A) in the state space H, we obtain (A2) and complete the proof 
of Theorem 2.4. 

Using Theorem 2.3 and Theorem 2.4, we can give the exponential stability of 
(2.1)-(2.5). 


Theorem 2.5. Let Y(0) € X be defined by (2.6) and set YP := PhY (0) be the initial 
data of (2.12)-(2.17), then Y} convergent to the initial value Y (0) of (2.1)-(2.5) in 
the sense of 


EYP + Y(0), as h + 0, 
and the semigroup T(t) is exponentially stable, i.e., 
ITŒY (O)\lx < Me™ |Y (0) |x, Y Y(0) €X, t20, 


for some constant M' and w is the same as in Theorem 2.3. 


Proof. As h — 0, EnY;? convergent to Y (0) since En and P, satisfy (A2). For any 
€ > 0, (b) of Theorem 2.4 implies that there exists he € (0,1) such that 


||EnTh (t) PRY (0) — THY Olly < € 
for all h < hz. However, it follows from (2.28) and (A1) that 
En Th (t) PRY (O)llx < M'e [Y (0)llx; for h < he, 
with M’ = M Mı M2. Therefore, combining the above two inequalities, we obtain 


T(L)Y 0)lly < Me™ +e. 


By letting € — 0, we derive that T(t) is exponentially stable. 
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3. The hybrid system with nonlinear feedback 

In this section, we cope with the hybrid system (1.1) and (1.7)-(1.8) 
Onu(x,t) = Opzu(x,t) + alxjut(x, t), 
mO,,u(0, t) — 0,u(0, t) = F(u:(0,t)), 
ue (1, t) = —g(ue(1, t)), 


Using the auxiliary functions v(#,t) = u(x,t) and w(z,t) = uz(x,t), the above 
hybrid system can be rewritten as 


Oulx, t) = 0,w(ax,t) — a(x)v(t, x), (3.1) 
Opw(a,t) = zvl, t), (3.2) 
mo,v(0,t) = F(v(0,t)) + w(0,t), (3.3) 
w(1,t) = —g(v(1, 4), (3.4) 
v(z,0) = v°(x), w(x,0) = w(x), (3.5) 


in which v°(x) and w(x) are given in the last section. a(x), g(s) and F(s) are scalar 


functions satisfying the following hypothesises?°: 


(Hı) a(x) € C10, 1] and |a(zx)| < a1 with ay € (0,1) and ma, < 4; 

(H2) g(s) is continuous, nondecreasing, and g(0) = 0; 

(H3) F(s) is globally g-Lipschitz continuous, and F'(0) = 0; 

H4) a2\|s| < |g(s)| < ag|s|, for all s € C and some positive real numbers 
Q2 and a3; 

e (H5) The inequality a2/(1 + a3) — q > aı(4 — may) holds. 


We refer the reader to Ref. 26 for well-posedness properties of the closed-loop system 
(3.1)-(3.5). 


3.1. UES of (3.1)-(3.5) 


Under the above assumptions, we analyze the UES of (3.1)-(3.5). For this purpose, 
we use the average central-difference scheme of last section to discretize the spacial 
derivative of (3.1)-(3.5) and obtain 


vya (t) = ôswjy 3 (t) — ajv; 42 (t), (3.6) 
w, a(t) = be0;43 (t), (3.7) 
mug(t) = F(vo(t)) + wo(t), (3.8) 
wynt) = —g(un4ilt)), (3.9) 
vp (0) = vg, w(0) = we, (3.10) 


where v;(t), w;(t), ak, v} and w? have the same meaning as in the last section. For 
simplicity, we drop for readability the coordinate t of v;(t) and w,(t) from now on. 
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The energy functional, which is denoted by Ena (t), of (3.6)-(3.10) is defined by 
N 


h 
Baalt) = 59 [lon 


j=0 


2 
oy oe 


2 m 
|+ Fir, (3.11) 


which satisfies similar result as proposition 2.1. 
Proposition 3.1. The discrete energy functional E,,;,(t) satisfies 


malt) =—h X. aluj? — unsig(wn41) + voF (v0). (3.12) 
jhew 


Proof. Multiplying (3.6), (3.7), and (3.8) by hvjya, hwy42, and vo respectively 
and summing for j, we derive 


N 
a(t) =h 5 a|vj42/ + hX lo; błôzWjyi + Wz 41020542]. (3.13) 
jhew j=0 
However, by using (2.23), we have 
N 
hY oj4 250s 44 + Wj 42600442] = UNGWN 41 — VOWD. (3.14) 
j=0 


Thus, (3.8)-(3.9) and (3.13)-(3.14) implies that (3.12) holds. 


Remark 3.1. Because no prior knowledge on the signs of the functions a and F, 
the energy functional E,,(t) is no longer non-increasing along the trajectories of 
the system. This means that F(s) and a(x) can be model destabilizing boundary 
and interior anti-damping phenomenons. 


To present main result of this section, we introduce two useful lemmas. The first 
one comes from Lemma 2.2 of Ref. 13 or Lemma 2.1 of Ref. 30. The second one is 
discrete counterpart of (57) of Ref. 26 or (A3) of Ref. 5. 


Lemma 3.1. Let {ui} h", {vi}M' and {wi}Mh! be sequences consisting of real 


numbers, then we have 


N N 
1 1 
a So (wis — u) (vipi Wieser + wi) + ri Cr — u) O41 =U) (Ua — By) 
i=0 i=0 
ee ies 
F 1 X (uia + ui) (vipi — vi) (wii + Wi) + 1 X (uima + ui) (vipi + Vi) (Wid — wi) 


i=0 i=0 
= UN+1VN+1WN+1 — UoVowo 
Lemma 3.2. Let p(x) be affine, positive, and increasing, p' = p' (x), p; = p(x;), 
and L(t) be the Lyapunov functional which is defined by 


N 


En(t) = Earl) + hd | Pye 304 RM 43> (3.15) 
j=0 
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then for any T > 0 we have 


T 
Ln(t)lo 
T T T 
z0 —20n) f Enn (t)dt — | opine ina TE | bral Faraldi 
0 0 0 


T T 
+ | voF(vo)dt — (290 + Mma /2 + mazp'/2) | vo |" dt. (3.16) 
0 0 


Proof. It is easy to know 


T T 
ma 
Enr Ol < | a1 Enn(t)dt J 2 vol? t un+19(UN +1) — vo F (uo) dt (3.17) 
0 0 


holds by (3.11)-(3.12) and (Hı). Just as in the proof of Proposition 3.1, we multiply 
(3.6) and (3.7) by hpj,1wj,1 and hp;,1v;,1 and obtain 


T 
N 
hY Pj+}Vj+}Wj44 
j=0 r 
T T N 
=- f h J ayuj yj gat | h J Pj+h [rig g berg + Fat pwt dt. 
0 jhew 0 5=0 


(3.18) 


Applying Lemma 3.1 for the second term in the right-hand side of (3.18), we obtain 
N 
hY Prag lrjrg Serj 43 + Fn; 43%; 43] 
j=0 
nA 
a So be 542lloj4a O + wga] + 2onsallonaal? + ww sil?) 
j=0 


h3 Š 
= 2pollvol? + jwol?] = Yo ap; 43ll6erj 44? + [Semj43P] (319) 
j=0 


Then, combining (3.18) and (3.19) yields 


N E T 
hÝ Pitit Wja] S -f (p' — a4) Enn(t)dt 
j=0 o 2 
T T 
+2ønsa | lewsal? + bwxrsal?dt ~ (2p + mayp'/2) | lvol2dt, (3.20) 
0 0 


where Cauchy-Schwartz inequality is used for the first term of the right-hand side 
of (3.18), 620;41 =p! > 0, (3.11), and (H1) are applied in (3.19). At last, plugging 
(3.20) into (3.17), we derive (3.16) and finish the proof of this lemma. 
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With the aid of Lemma 3.2, we can derive UES of (3.6)-(3.10) under the hy- 
pothesises (Hı)-(ĦH5). 


Theorem 3.1. Assume that (H; )-(H;) hold, the discrete systems (3.6)-(3.10) are 
uniformly exponentially stable with respect to the energy Enn (t), i.e., there exist two 
positive constants M and w, (independent if h and solutions) such that 


Enn(t) < Me~“"*E,1,(0), for allt > 0. (3.21) 


Proof. Firstly, from (3.15) we know 
(1 — a1) Enn(t) < Ena(t) < (1 + 01) Enn(t) (3.22) 
since 
N 
1 1 1 
Jn X 82 pj5?0;52 w| < ar Ena(t). 
j=0 


Secondly, (H3) implies sF (s) < q|s|? for all s € R and 


T T T 
| voF (vo)dt — (200 + may /2 + mayp'/2) f |uo|?dt < =ma,/2 | |uo|7 dt, 
0 0 0 


(3.23) 
where 299 + may;p'/2 = q is assumed for a moment. By using (3.9) and (H4), we 
see that 

T T 
ons | lunyil? + [wnt dt < 2PN41(1 + 03) f lun+1|?dt. (3.24) 
0 0 


Moreover, combing (H2) and (H4) yields 


T T 
-f vy+1g(vn+1)dt < -o> | |un4i|?dt. (3.25) 
0 0 


(S) Finally, assuming py+1 = a2/2(1 + a3), we can know that the affine function 
p(x) is uniquely defined by noticing py+1 = po + p’ and collecting above relations 
on p 


2po + marp'/2 =q, 

a2 
2(1 + a3)’ 
po and p’ are uniquely existed since ma, < 4 and p(x) is unique defined by po and 
p’ and indeed increasing with p’ = 2[a2/(1+0a%)—q|/(4—maz) > 0 (see H5). Now, 
combing (3.16) and (3.22)-(3.25), we obtain 


pot p= 


T 
Lant) < -w f Lnn(t)dt, for all T > 0. 
0 


where 
2la2/(1 + 03) — q]/(4 — may) — 201 


w= 
l-a, 
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By applying (Hı), (H5), and Grönwall lemma, we obtain the UES of (3.21). The 
proof of the Theorem 3.1 is complete. 


Remark 3.2. We can perform convergence analysis for the discrete system (3.6)- 
(3.10) by using the methods of Ref. 13, Ref. 5 or Ref. 24. But we delete it since the 
methods are very standard and the dynamical boundary conditions and nonlinear 
feedback do not effect the convergence. Moreover, we know that the discrete energy 
functional Enn also is convergent to the energy functional of (3.1)-(3.5) by similar 
result Theorem 5.1 of Ref. 13. Therefore, the nonlinear counterpart of Theorem 2.5 
is still right by the same analysis. That is to say (3.1)-(3.5) is exponentially stable 
T under the hypothesis of (H1)-(H5). 


3.2. Numerical simulations 


Firstly, we rearrange (3.6)-(3.10) into vectorial formulation. The unknown of (3.6)- 
(3.10) is Yna(t) = (vn (t), un +i (t), wa(t))*. up(t), wrt), Bn4i, MN+ı1, El; and Na 
are the same as those of subsection 2.2. Moreover, set a(x) = 1/10 cos(mx), q = 1/10, 


a2 = 1, a3 = 2, ey41 = (0,0,--- ,0,1) E CWD to be row vector. 
Secondly, in light of hypothesises (H2)-(Hs), for s € C we choose nonlinear 
functions 


1 
F(s) = or Vat 6: s>0, 
-7 + VI + 98, s <0, 


and 
© -14+2s+—4,s>0 
= _ st1?" = ™ 2 
= g(s) 1428-5, s<0. (3-26) 


Q F(s) clearly satisfies (H3) due to 


d 1 
—F(s)=q- ————__,, <4, (3.27) 
ds (ls +712) 


which implies that F(s) is globally q-Lipschitz continuous. g(s) is obviously contin- 
uous and g(0) = 0. Moreover, g(s) satisfies (H2) and (H4) since 


d 1 


— qls) = 2- ——__... 3.28 
49) =2- Tap (3.28) 


At last, we choose m = 39 and it is easy to verify that a; = 1/10, q = 1/10, ag = 1, 
a3 = 2 and m satisfy (H5). Thus, using the mean value Theorem, F(0) = 0, and 
g(0) = 0, we replace F(vo(t)) and g(un4i(t)) by 


7 F(vo(t)) and $ glon), 
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respectively. Therefore, we obtain the quasi-linearized system of (3.6)-(3.9) 


vya (t) = d2wj4a(t) — ajvj (t+ - (3.29) 
w (t) = Seva (t); (3.30) 
h(t) = EF (volt) )volt) + wolt), (3.31) 
wnli) = -Eglon Oenslt). (3.32) 


In order to perform numerical simulations and give ideal approximating of 


< F(vo(t)) and <alowsr(t)); 


we replace s by the positive random number ranging from 0 to 10000, which is 
denoted by r, in (3.27) and (3.28). We define the matrixes Bnn and Mnn by 


Bn eN41 0 
Bat) 2 0 0 
2 1 * 
0 (mi) engi BN+1 
and 
i 2 BniiNa ( hans 2+ — j| eN41 “MN 
= h 

Mnh = h ghey paye 0 0 
—Myy1 eN+1 0 


Because Bnn is evidently invertible, the linearized system (3.29)-(3.32) are then 
obtained 


(S) BanVi p(t) = MahYnn (t) @ Ya) = Anh Yn (t), with Ann := B7} Mnn. (8.33) 


Similarly as in subsection 2.2, one replace Bnn by identity operator and the lin- 
earized classical finite difference scheme of (3.1)-(3.4) is easily restored from (3.33), 
i.e., 


Ynh (t) = MnhYnh (t). 


Here we explain the significance of the discrete scheme (3.33) by plotting two 
kinds of figures. Figure 3 gives the maximal real parts of Ann and Mnn for N = 40: 
5 : 400. Figure 4 depicts the distributions of the eigenvalues of Mnn with N = 500. 
We see that the real parts of the eigenvalues of Mnn approach to zero form above 
and those of Anp approach to a negative number from these figures. Numerical 
simulation results show that the classical finite difference scheme of (3.1)-(3.4) is 
not uniformly exponentially stable and even not stable. This conclusion is consistent 
with it’s non-dissipativity or anti-damping. 
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4. Concluding remarks 


The exponential stabilities of two hybrid systems are investigated from different 
viewpoint as before. More precisely, semi-discretized systems and their UESs are 
first derived. The exponential stability of the hybrid system follows from the conver- 
gence of numerical approximating scheme and the result of UES. Looking back at 
the Subsection 2.2, one may find out that the linear feedback —au,(0,t) + &aust(0, t) 
play key role for UES of discrete hybrid system with linear feedback. However, the 
local viscosity term a(x)uz(#,t) is trivial to the main result. This phenomenon also 
appears in the hybrid system with nonlinear feedback of Section 3. This manifest 
that velocity stabilizer or velocity plugging angular velocity stabilizer is power- 
ful. We adopt two approaches, namely the frequency domain method and discrete 
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multiplier method, to verify the UES of corresponding system. In the future, it is 
interesting to apply frequency domain verification method developed in this paper 
to other systems with boundary control. Moreover, we only consider one case of 
the nonlinear hybrid system of Ref. 26, it is worth investigating other cases in the 
further research. 


Acknowledgment 


The authors would like to thank the anonymous reviewers for their very careful 
reading and valuable comments and suggestions to improve the manuscript. 


References 


1. 


10. 


11. 


12. 


13. 


F. Abdallah, et al., Uniformly exponentially or polynomially stable approximations 
for second order evolution equations and some applications, ESAIM: COCV. 19 (2) 
(2013) 844-887. https://doi-org/10.1051/cocv/2012036 

H.T. Banks, K. Ito, C. Wang, Exponentially stable approximations of weakly damped 
wave equations, in Estimation and control of distributed parameter systems (Basel: 
Birkhauser, 1991). 

P. Freitas, E. Zuazua, Stability results for the wave equation with indefinite damping, 
Journal of Differential Equation 132 (2) (1996) 338-353. 

https: //doi.org/10.1006/jdeq.1996.0183 

B.Z. Guo, F.F. Jin, Arbitrary decay rate for two connected strings with joint anti- 
damping by boundary output feedback, Automatica 46 (7) (2010) 1203-1209. 
https://doi.org/10.1016/j.automatica.2010.03.019 

B.Z. Guo, B.B. Xu, A semi-discrete finite difference method to uniform stabilization 
of wave equation with local viscosity, IFAC Journal of Systems and Control 13 (2020) 
101000. https: //doi.org/10.1016/j.ifacsc.2020.100100 

B.Z. Guo, C.Z. Xu, On the spectrum-determined growth condition of a vibration cable 
with a tip mass, IEEE Transactions on Automatic Control 45 (1) (2000) 89-93. 
https://doi.org/10.1109/9.827360 

B.Z. Guo, F. Zheng, Frequency energy multiplier approach to uniform exponential 
stability analysis of semi-discrete scheme for a Schrödinger equation under boundary 
feedback, Preprint, Hainan University, 2022. 

F. Hassine, Rapid exponential stabilization of a 1-d transmission wave equation with 
in-domain anti-damping, Asian Journal of Control 19 (6) (2017) 2017-2027. 
https://doi.org/10.1002/asjc.1509 

J.A. Infante, E. Zuazua, Boundary observability for the space semi-discretizations of 
the 1-d wave equation, Mathematical Modelling and Numerical Analysis, 33 (2) (1999) 
407-438. https://doi-org/10.1051/m2an:1999123 

K. Ito, F. Kappel, The Trotter-Kato Theorem and Approximation of PDEs, Math- 
ematics of Computation 221 (8) (1998) 21-44. https://doi.org/10.1090/S0025-5718- 
98-00915-6 

B. Jacob, H. Zwart, Linear Port-Hamiltonian System on Infinite-dimensional Space 
(Basel, Switzerland: Springer-Verlag, 2012). 

R. Kress, Numerical Analysis, Graduate Texts in Mathematics 181 (New York: 
Springer-Verlag, 1998). 

J.K. Liu, B.Z. Guo, A New Semi-discretized Order Reduction Finite Difference Scheme 
for Uniform Approximation of 1-D Wave Equation, SIAM J. Control Optim. 58 (4) 
(2020) 2256-2287. https://doi.org/10.1137/19m1246535 


January 16, 2023 8:23 WSPC/INSTRUCTION FILE ws-m3as 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22: 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


30. 


Uniform exponential stability approximations for two hybrid systems 25 


J.K. Liu, B.Z. Guo, A novel semi-discrete scheme preserving uniformly exponential 
stability for an Euler-Bernoulli beam, Systems & Control Letters 160 (2019) 104518. 
https://doi.org/10.1016/j.sysconle.2019.104518 

J.K. Liu, R.Q. Hao, B.Z. Guo, Order reduction-based uniform approximation of expo- 
nential stability for one-dimensional Schrödinger equation, Systems & Control Letters 
160 (2022) 105136. https: //doi.org/10.1016/j.sysconle.2022.105136 

Z.Y. Liu, S.M. Zheng, Uniform Exponential Stability and Approximation in Control 
of a Thermoelastic System, SIAM J. Control Optim. 32 (5) (1994) 1226-1246. 
https: //doi.org/10.1137/S0363012991219006 

Z. Liu, S. Zheng, Semigroups associated with dissipative systems, 398 of Chapman 
& Hall/CRC Research Notes in Mathematics (Chapman & Hall/CRC, Boca Raton 
1999). 

A. Macchelli, Y.L. Gorrec, H. Ramirez, Exponential stabilization of port-Hamiltonian 
boundary control systems via energy shaping, IEEE Transactions on automatic control 
65 (10) (2020) 4440-4447. https: //doi.org/10.1109/TAC.2020.3004798 

A. Macchelli, Y.L. Gorrec, Y. Wu, H. Ramirez, Energy-based control of a wave equa- 
tion with boundary anti-damping, [FAC-PapersOnLine 53 (2) (2020) 7740-7745. 
https://doi.org/10.1016/j.ifacol.2020.12.1527 

Z.D. Mei, Output feedback exponential stabilization for a 1-d wave PDE with dynamic 
boundary, J. Math. Anal. Appl. 508 (1) (2022) 125860. 
https://doi.org/10.1016/j.jmaa.2021.125860 

O. Morgul, B.P. Rao, F. Conrad, On the stabilization of a cable with a tip mass, JEEE 
Transactions on automatic control 39 (10) (1994) 2140-2145. 
https://doi.org/10.1109/9.328811 

A. Smyshlyaev, M. Krstic, Boundary control of an anti-stable wave equation with 
anti-damping on the uncontrolled boundary, Systems & Control Letters 58 (8) (2009) 
617-623. https://doi.org/10.1016/j.sysconle.2009.04.005 

L.T. Tebou, E. ZuaZua, Uniform exponential long time decay for the space semi- 
discretization of a locally damped wave equation via an artificial numerical viscosity, 
Numerische Mathematik 95 (4) (2003) 563-598. https://doi.org/10.1007/s00211-002- 
0442-9 

L.T. Tebou, E. ZuaZua, Uniform boundary stabilization of the finite difference space 
discretization of the 1-d wave equation, Advances in Computational Mathematics 26 
(1) (2007) 337-365. https: //doi.org/10.1007/s10444-004-7629-9. 

A. Terrand-Jeanne, et al., Regulation of inhomogeneous drilling model with a P-I 
controller, IEEE Transactions on Automatic Control, 65 (1) (2020) 58-71. 
https://doi.org/10.1109/TAC.2019.2907792. 

N. Vanspranghe, F. Ferrante, C. Prieur, Velocity stabilization of a wave equation with 
a nonlinear dynamic boundary condition, IEEE Transactions on Automatic Control, 
to be published. 

X.F. Wang, W.L. Xue, Y. He, F. Zheng, Uniformly exponentially stable approximations 
for Timoshenko beams, Submitted. 

Y.L. Zhang, M. Zhu, D.H. Li, J.M. Wang, Stabilization of two coupled wave equa- 
tions with joint anti-damping and non-collocated control, Automatica 135 (1) (2022) 
109995. https: //doi.org/10.1016/j.automatica.2021.109995 

F. Zheng, S.J. Zhang, H.K. Wang, B.Z. Guo, The exponential stabilization of a heat- 
wave coupled system and its approximation. J.Math. Anal. Appl. 521 (2023) 126927. 
https: //doi.org/10.1016/j.jmaa.2022.126927 

F. Zheng, H. Zhou, State reconstruction of the wave equation with general viscosity 
and non-collocated observation and control, J. Math. Anal. 502 (2020) 125257. 


202301.00191v1 


chinaXiv 


January 16, 2023 8:23 WSPC/INSTRUCTION FILE ws-ms3as 


26 F. Zheng, L. Zhang, Y. He 


https://doi.org/10.1016/j.jmaa.2021.125257 

31. E. Zuazua, Propagation, observation, and control of waves approximated by finite 
difference methods, SIAM Rev. 47 (2005) 197-243. 
https: //doi.org/10.1137 /S0036144503432862 


